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Abstract 

We carry out a detailed analysis of the nonminimal supersymmetric standard model with lepton 
number violation. The model contains a unique trilinear lepton number violating term in the 
superpotential which can give rise to neutrino masses at the tree level. We search for the gauged 
discrete symmetries realized by cyclic groups Zjy which preserve the structure of the associated 
trilinear superpotential of this model, and which satisfy the constraints of the anomaly cancellation. 
The implications of this trilinear lepton number violating term in the superpotential and the 
associated soft supersymmetry breaking term on the phenomenology of the light neutrino masses 
and mixing is studied in detail. We evaluate the tree and loop level contributions to the neutrino 
mass matrix in this model. We search for possible suppression mechanism which could explain 
large hierarchies and maximal mixing angles. 
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I. INTRODUCTION 



Supersymmetry [1] is at present the only known framework in which the Higgs sector of 
the Standard Model (SM), so crucial for its internal consistency is natural. A much favored 
implementation of the idea of supersymmetry at low energies is the minimal supersymmetric 
standard model (MSSM), which is obtained by doubling the number of states of SM, and 
introducing a second Higgs doublet (with opposite hypercharge to the Standard Model 
Higgs doublet) to generate masses for all the SM fermions and to cancel the triangle gauge 
anomalies. However, the MSSM suffers from the so-called /i problem associated with the 
bilinear term connecting the two Higgs doublet superfields H u and H d in the superpotential. 
An elegant solution to this problem is to postulate the existence of a chiral electroweak 
gauge singlet superfield S, and couple it to the two Higgs doublet superfields H u and Hd via 
a dimensionless trilinear term XH u HdS in the superpotential. When the scalar component of 
the singlet superfield S obtains a vacuum expectation value, a bilinear term \H u H d < S > 
involving the two Higgs doublets is naturally generated. Furthermore, when this scalar 
component of the chiral singlet superfield S acquires a vacuum expectation value of the order 
of the SU (2) L xU (l)y breaking scale, it gives rise to an effective value of /i, /j, e ff = A < S >, 
of the order of electroweak scale. However, the inclusion of the singlet superfield leads to 
additional trilinear superpotential coupling («;/3)S' 3 in the model, the so called nonminimal, 
or next-to-minimal [2-6], supersymmetric standard model (NMSSM). The absence of H u Hd 
term, and the absence of tadpole and mass couplings, S and S 2 in the NMSSM is made 
natural by postulating a suitable discrete symmetry. The NMSSM is attractive on account 
of the simple resolution it offers to the /i-problem, and of the scale invariance of its classical 
action in the supersymmetric limit. Since no dimensional supersymmetric parameters are 
present in the superpotential of NMSSM, it is the simplest supersymmetric extension of the 
Standard Model in which the electroweak scale originates from the supersymmetry breaking 
scale only. Its enlarged Higgs sector may help in relaxing the fine-tuning and little hierarchy 
problems of the MSSM [7], thereby opening new perspectives for the Higgs boson searches 
at high energy colliders [8-11], and for dark matter searches [12]. 

Since supersymmetry requires the introduction of superpartners of all known particles in 
the SM, which transform in an identical manner under the gauge group, there are additional 
Yukawa couplings in supersymmetric models which violate [13] baryon number (B) or lepton 
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number (L). In the minimal supersymmetric standard model there are both bilinear and 
trilinear lepton number violating Yukawa terms in the superpotential. There are also trilinear 
baryon number violating Yukawa terms in the superpotential. All these terms are allowed by 
renormalizability, gauge invariance and supersymmetry. In MSSM, a discrete symmetry [14] 
called i?-parity (Ftp) is invoked to eliminate these B and L violating Yukawa couplings. 
However, the assumption of i?-parity conservation at the level of low energy supersymmetry 
appears to be ad hoc, since it is not required for the internal consistency of supersymmetric 
models. 

If we do not postulate R p conservation, then there are baryon and lepton number number 
violating terms in the superpotential of NMSSM as well. What is perhaps interesting is the 
presence of an additional lepton number violating trilinear superpotential coupling [15, 16] 
in this model which has no analog in the MSSM with baryon and lepton number violation. It 
is, therefore, important to study the implications of this additional lepton number violating 
trilinear interaction term in the superpotential of NMSSM, contrast the situation with MSSM 
with lepton number violation, and pin down the possible differences with its predictions. 

One of the far reaching implications of the lepton number violating couplings in NMSSM 
concerns the physics of light neutrino states. In identifying the dominant contributions to 
the neutrino masses, and suppression mechanisms, one must compare with the situation 
that obtains in MSSM with bilinear lepton number violation. In NMSSM, the three light 
neutrinos mix with SU(2) L x U(l)y nonsinglet gaugino and Higgsino fields as well as the 
gauge singlet fermionic component of S, the singlino (S). The resulting 8x8 mass matrix 
of the neutrino-gaugino-Higgsino-singlino has a see-saw structure, which leads to a sepa- 
rable rank one effective mass matrix for the light neutrinos, implying the presence of a 
single massive Majorana neutrino. At one-loop order, there occur two main mechanisms 
for generating masses for the Majorana neutrinos. One of these involves only the matter 
interactions. The second mechanism involves matter interactions in combination with the 
gauge interactions and propagation of neutralinos and mixed sneutrino-Higgs boson sys- 
tem, whose contribution depends sensitively on the soft supersymmetry breaking couplings. 
While both these mechanisms have the ability to generate masses for the Majorana neutri- 
nos, the latter one, initially proposed in the context of MSSM by Grossman and Haber [17], 
is expected to dominate. In the case of MSSM with bilinear lepton number (or R p ) violation, 
the tree and one-loop contributions to the neutrino masses, and their ability to reproduce 
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the experimental observations, have been extensively discussed in the literature [18-26]. 

In this paper we carry out a detailed investigation of the nonminimal supersymmetric 
standard model with lepton number violation. Since the NMSSM has a unique lepton 
number violating trilinear coupling term in the superpotential, one of the issues we want to 
address concerns the implications of the neutrino masses and mixing for the NMSSM with 
such a lepton number violating term. Our purpose is to pin down the features specific to this 
version of the NMSSM and extract constraints implied by the comparison with experimental 
data. We compare and contrast the situation in NMSSM with that of MSSM with bilinear 
i?-parity violation (RPV) [24]. Despite the presence of the singlino, and its mixing with 
the neutrinos, a light mass Majorana neutrino appears at the tree level. This is due to 
the constrained nature of the couplings in the model. Nevertheless, as in the MSSM with 
lepton number violation, one-loop contributions play an important role in determining the 
neutrino mass spectrum. The ability to reproduce experimental observations is expected to 
set useful constraints on the Higgs sector parameters of the NMSSM. The situation differs 
from the one that arises in the see-saw mechanism or the bilinear lepton number violation in 
MSSM in that no dimensional mass parameters (large or small) are introduced. The neutrino 
Majorana masses arise from dimensionless Yukawa couplings. However, despite the presence 
of a gauge singlet fermion that could play the role of a sterile neutrino, whether an ultra light 
singlino mode, compatible with the the cosmological bound on the summed mass of light 
neutrinos, J2u m u < 10 eV, does indeed occur is at variance with the physical constraints on 
the NMSSM which rule out the possibility that the lightest mode in the massive neutralino 
sector lies below O(50) GeV. Thus, the understanding of neutrino physics provided by the 
NMSSM with lepton number violation contrasts with that proposed in models using the 
compactification moduli superfields [27] or axion superfields [28] coupled gravitationally to 
the observable sector modes. 

This paper is organized as follows. We begin in Section II with a discussion of the 
general structure of the superpotential and soft supersymmetry breaking interactions in the 
nonminimal supersymmetric standard model with lepton number violation. In this Section 
we also discuss the local Z N cyclic symmetries which can protect the NMSSM against B or 
L, or combined B and L, number violating superpotential couplings (Subsection II B). We 
further elaborate on the general approach to analyze the gauged cyclic group symmetries 
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in the Appendix A to the paper. In Section III we derive the tree level light neutrino 
mass spectrum that arises in this model (subsection III A), and then obtain the one-loop 
radiative corrections to the mass spectrum in subsection III B. In Section IV we present a 
general discussion of the predictions from this model, which are based on the consideration 
of Abelian horizontal symmetries for the flavor structure of effective couplings. Finally in 
Section V we summarize our results and conclusions. 

II. NMSSM WITH BARYON AND LEPTON NUMBER VIOLATION 
A. The superpotential 

In this section we recall the basic features of NMSSM with baryon and lepton number 
violation, and establish our notations and conventions. The superpotential of the model is 
written as 

Wnmssm = {hu)abQVJ rH u + {hD)abQ a LD R H d + (fl e) abL L E R H d + \SH d H u — -S* 3 , (II. 1) 

where L, Q, E, D, U denote the lepton and quark doublets, and anti-lepton singlet, d-type 
anti-quark singlet and u-type anti-quark singlet, respectively. In Eq. (II. 1), (hjj) ab , {h D ) ab 
and {hE)ab are the Yukawa coupling matrices, with a, b, c as the generation indices. Gauge 
invariance, supersymmetry and renormalizability allow the addition of the following L and 
B violating terms to the superpotential (II. 1): 

W L = \ a L a H u S + l -\ ahc L a L L L E R + \' abc L a L Q b L D R , (11.2) 

W B = \KbcD R D b R U c w (II.3) 

where the notation [15, 16] is standard. We note that there is an additional L-violating 
term with the dimensionless Yukawa coupling A a in (II. 2) which does not have an analogue 
in the MSSM. This term can be rotated away into the R-parity conserving term \SH u H d 
via an SU(4) rotation between the superfields H d and L a . However, this rotation must be 
performed at some energy scale, and the term is regenerated through the renormalization 
group equations. The Yukawa couplings \ a b c and \" bc are antisymmetric in their first two 
indices due to SU(2) L and SU(3)c group symmetries, respectively. 

The supersymmetric part of the Lagrangian of NMSSM with baryon and lepton number 
violation can be obtained from the superpotential (II. 1), (II. 2) and (II. 3) by the standard 
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procedure. In addition to this supersymmetric Lagrangian, there are soft supersymmetry 
breaking terms which include soft masses for all scalars, gaugino masses, and trilinear scalar 
couplings, respectively. These can be written as 



soft 



-c 



soft 



M Q 2 Q a *Q h L + M u 2 U R U R + M$ 2 D R D R + M2 b2 L a L *L b L + M a E b2 E R E R 
+ m 2 H H* d H d + m 2 Hu H* u H u + m 2 s S*S 
+ ^M s \ s \ s + ^M 2 \ w \ w + ^MxA'A' 

„ ^_b _ ^_b 

+ {Au) ah (hu) ab Q a L U R H U + (A D ) ab {h D ) ab Q a L D R H d 
+ {A E ) ab {h E ) ab L a L E R H d - A x XSH d H u - ^kS 3 



+ 



~ ~ , ~ c 



~( A \ a )^a,L a L H u S + -(A X )abcKbcL a L L L E R + (A X ')abcKbc L LQL D R 



+ 



~ a ~ ft ~ c 
(^A")abcA" 6c -D i?J D R f/ jR 



+ H. c, 



(II.4) 



where a tilde over a matter chiral superfield denotes its scalar component, and the notation 
for the scalar component of the Higgs superfield is the same as that of the corresponding 
superfield. We note that the soft supersymmetry breaking gaugino masses have been de- 
noted by Mi,M 2 , and M s corresponding to the gauge groups U(l)y, SU(2) L , and SU(3)c, 
respectively. We have chosen the sign conventions for the soft trilinear couplings involving 
the gauge singlet field in Eq.(II.4) which are different from those used in Ref [16]. 

The dimension-4 terms in the superpotentials (II. 2) and (II. 3) are the most dangerous 
terms for nucleon decay, and some of them must be suppressed. This leads to constraints [29] 
on the different couplings \ abc , \' abc , and \" bc , but considerable freedom remains for the 
various B and L violating couplings. Furthermore, there are dimension-5 operators [16] 
which may lead to nucleon decay suppressed by 1/M, where M is some large mass scale 
at which the B and L violation beyond that of NMSSM (and MSSM) comes into play. 
Some of these dimension-5 operators may also lead to unacceptable nucleon decay if their 
coefficients are of order unity, and therefore must be suppressed. We shall not consider 
here the dimension-5 operators, but instead concentrate on the dimension-4 lepton number 
violating terms (II. 2) only. 

As noted above, there is an additional trilinear L-violating term with the dimensionless 
Yukawa coupling A a in (II. 2) which does not have an analogue in the MSSM. This term can 
be rotated away into the R-parity conserving term \SH u H d via an SU (4) rotation between 
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the superfields Hd and L a . However, this rotation must be performed at some energy scale, 
and the term is regenerated through the renormalization group equation [15, 16] 



<1Xa ' '(:*/,? • lr • 4A 2 • -2,r • 1A 2 • A 2 33 • SA^) A 3 



din// IQtt 2 

+3h b X\' 333 - (Sgj + -gjfj A 3 



(II.5) 



where (^17)33 ~ h t , (h D ) 33 ~ h b , (h L ) 33 ~ /i r . Here £/ 2 , #1 are the SU(2) L and £/(l)y gauge 
couplings, and the meaning of other quantities is obvious. For simplicity, we have retained 
only the highest generation trilinear couplings in (II. 5). 

It is important to point out two distinctive features of the present model relative to the 
MSSM. First, in NMSSM no distinction is made between the bilinear and trilinear lepton 
number violation, since the bilinear Lagrangian terms, fiHdH u + /iiLiH u , [fi = A < S > 
, Hi — Aj < S >], can arise as effective couplings once the singlet scalar field component of 
S acquires a finite VEV. The wide hierarchies between the lepton number conserving and 
violating couplings, of the expected size \/\ = fiij fi ~ 1CT 6 , arise from the hierarchies of 
the dimensionless Yukawa couplings. Second, the lepton number violating trilinear operator 
LiH u S has the ability to radiatively induce other trilinear lepton number violating couplings, 
which is precluded for the bilinear operator LiH u . This property may be used to justify a 
scheme where naturally suppressed trilinear couplings A^, \'^ k occur as a result of being set 
to zero at some large mass scale (gauge unification scale), and receive small finite radiative 
corrections from the gauge singlet coupling XiLiH d S. This possibility can be established 
on a quantitative basis by examining the one-loop renormalization group equations for the 
trilinear coupling with maximal number of third generation indices [16] 

^f^Q = A 2 33[4^ + 4AL + 3A; 2 33 + A^-(^ 1 2 + 3, 2 2 )], 

r)\' 

(47r)2 ^o?Q = A ' 333 ^ + 6hl + K + ~ Xl + AL + 6X ™ 
16 7 
-(y<?3 + 3<? 2 2 + -<? 1 2 )]+MA 3 . (116) 

Looking for solutions of these equations with suitably large values of the trilinear couplings 
at some large mass scale, for instance, gauge coupling unification scale, might reveal the 
presence of infrared fixed points which would then serve as upper bounds on the weak scale 
values of these couplings. 
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To establish the lepton number violating nonminimal supersymmetric standard model 
on a firmer basis, it is important to determine whether there are discrete symmetries re- 
specting the postulated interaction superpotential which can be regarded as local or gauged 
symmetries [30] obeying the anomaly cancellation conditions. As is known, the gauged dis- 
crete symmetries enjoy a natural protection against breaking by non-perturbative quantum 
effects initiated by the gravitational interactions, and against the emergence of massless 
Nambu-Goldstone bosons from the spontaneous symmetry breaking. Another advantage 
lies in evading the cosmological domain wall problem by the removal of classical domain 
wall solutions as a result of the gauge equivalence of degenerate vacua. We recall that if 
stable domain wall solutions were present, the production of cosmic domain walls in the 
early universe would result in a contribution to the present day mass density of the universe 
which exceeds the critical density. The case of Green-Schwarz (GS) anomalous gauged dis- 
crete symmetries is special in that although domain wall solutions do exist in gauge field 
theories satisfying such global type symmetries, the instanton tunneling effects present in 
these theories lift the degeneracy of vacua so as to render the solutions unstable [31, 32]. The 
above resolution applies independently of the familiar one invoking the domain wall dilution 
during a cosmic inflation era. We also note that the discrete gauge symmetries have been 
used in connection with various naturalness problems, such as the doublet-triplet splitting 
in unified gauge models, the stabilization of axion symmetries, or the construction of flavor 
symmetries realized via Froggatt-Nielsen mechanism [33, 34]. Following the work by Ibanez 
and Ross [35, 36], we are led here to consider the so-called generalized parities (GPs) of the 
NMSSM which forbid part or all of the dangerous couplings at the renormalizable level [37]. 
In the following Subsection we shall describe the construction of the generalized parities for 
the NMSSM in detail. 

B. Discrete symmetries 

In this Subsection we consider the Z N cyclic local symmetries which can protect the 
NMSSM against B, or L, or combined B and L number violating superpotential couplings. 
Demanding consistency with the anomaly cancellation conditions sets highly non-trivial 
constraints on the generalized baryon, lepton and matter parity symmetries (designated as 
GBP, GLP and GMP) for the appropriate superpotential. These are examined by making 
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use of the approach of Ibahez and Ross [35, 36], which is developed in Appendix A for the 
NMSSM. Before presenting our results, we shall outline the problem in general and introduce 
our notations. 

Besides the regular R-parity conserving (RPC) trilinear matter-Higgs Yukawa couplings, 
QH U U C , QH d D c , LH d E c , and the dangerous R-parity violating (RPV), and B and L 
violating, couplings U C D C D C and LLE C , LQD C , the renormalizable superpotential of the 
NMSSM includes the trilinear couplings H d H u S and LH U S, but excludes the dimensional 
superpotential couplings S, S 2 and H d H u . Thus, in addition to the weak hypercharge 
Y, the regular couplings conserve three U(l) charges. A convenient basis for the gen- 
erators of the corresponding 3-dimensional vector space is given by the U(l) charges 
R, A, L, which are defined in the table of Appendix A. The cyclic group elements 
R = e iaRR , A = e lOLA , and L = e %OLhL are defined by restricting the complex phase 
rotation angles to the values an = 2nm/N, aA = 2irn/N, cxl = 2np/N, with m,n,p 
being integers. The generators of the independent Z N multiplicative symmetry groups 
are thus of the form, g = R m A n L p = g n PQ R m ~ 2n 'If ', [g PQ = R 2 A] with specific (modulo 
N) relations linking the three integers m, n, p, depending on the allowed subset of dan- 
gerous couplings. For the ordinary symmetries (O), the charges are given explicitly by 
9(Q) = °> 9(U C ) = -m, g(D c ) = m- n, g(L) = -n - p, g(E c ) = m + p, g(H d ) = 
—m+n, g(H u ) = m, g(S) = mx+ny+pz. The charges x, y, z assigned to the gauge singlet S 
must obey the selection rules, x+y+z ^ 0, 2(x+y+z) ^ 0, 3(x+y+z) = 0, x+y+z+n = 0. 
A special role is played by the Peccei-Quinn like charge, (PQ) = 2R + A, which has a finite 
color group anomaly and is conserved by all the NMSSM couplings except S 3 , so a massless 
axion mode would arise only if the cubic coupling were absent, k — 0. The general form 
of the GBP, GLP and GMP generators is displayed in eq. (A. 2) of Appendix A. The soft 
supersymmetry breaking terms, which are generated via couplings to the goldstino spurion 
superfield X of form, V so f t ~ [AW] F + H. c. = F X W + H. c, are automatically protected 
for ordinary symmetries (O), and this protection remains valid for the R symmetries (R) as 
well, provided one assigns the R charge Q(X) = 0, and hence Q(Fx) = 2. 

Having classified the Abelian GP generators of fixed order N in terms of the three integers 
(m, n, p), our next task is to select the solutions satisfying the quantum anomaly conditions. 
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TABLE I: The solutions for the generalized flavor blind Zn parity symmetries of the NMSSM 
which cancel the mixed gauge anomalies ^3,2,1 only. We have made the choice N g = 3, N2h = 1 
for the number of quark and lepton generations and Higgs boson pairs, respectively, and k\ = | 
for the hypercharge normalization. The cancellation conditions for the gravitational and non-linear 
anomalies, A grav , Az?, Az» are not obeyed in general. The four rows give the solutions for the 
GBP, GLP, and GMP generators for symmetries of four distinct types: ordinary anomaly free 
(O), ordinary Green-Schwarz anomalous (O/GS), R symmetry anomaly free (R), and anomalous 
R symmetry (R/GS), respectively. The entries display integer (modulo N) parameters (m,n,p) 
for the Zn generator, and in the suffix we have given the values of the finite anomaly coefficients 
Agrav, Az2. For simplicity, we have not quoted the generally finite chiral anomaly coefficient A ^3. 
For the GMP symmetries, additional solutions of same order N as those quoted in the Table below 
also arise, as discussed in the text. 



Z N Type 


N 


GBP 


GLP 


GMP 


O 


9 


(l,3,4)_ 5 l,-87 
(5,6,2)_ 75 _34 8 
(2, 6, 8)_i02,-348 


(6,3,4)_36_ii 7 
(3, 6, 2)_ 8 i 5 _252 
(3, 6, 8)_99 5 _360 


(0,3,4)_ 54i _8i 
(2,3,4)_ 48i _ 93 
(3,3, 4)„ 45i _99 


O/GS 


7 


(0, 6, 2)_i26,-108 


(5, 6, 2)_i4i i _34 8 


(1, 6, 2)_i29,-l56, (2, 6, 2)_i32 -240 
(3, 6, 2)_i 35) _252 , (4, 6, 2)_i38 -300 


R 


12 


(8, 6, 6)_i3i 5 _560 


(2, 6, 6)_i49 i _320 


(0, 6, 6)_i55 _240 , (1,6, 6)_i52 -280 


R/GS 


11 


(5, 6, 2)_i 60i _368 


(3, 6, 2)_i54 5 _240 
1 n 


(0, 6, 2)_i45 -4 8 , (1, 6, 2)_i48 j _ii2 



These are expressed in terms of an over-determined system of linear and nonlinear equations 
for m, n, p. We have developed a numerical program to solve the anomaly cancellation 
conditions for the generalized parity generators (GBP, GLP, GMP) of four types. Unless 
stated otherwise, the search has been restricted to the case involving three quark and lepton 
generations, N g = 3, and a single pair of Higgs boson doublets, N 2 h = 1- While our 
presentation of results will be limited to cyclic groups of order iV < 15, we note that higher 
order solutions occur at integer multiples of the low order solutions, and as such they do not 
reveal novel features. 

If we demand that all the anomaly constraints are satisfied, then we find that, in general, 
no solutions exist, even if one is willing to push the search to high enough group orders, say 
N < 30. However, one must realize that the various anomaly cancellation conditions need 
not all be placed on the same footing. The linear anomalies A3, A2, Ai, A gra v have an 
obvious priority over the others, to the extent that these identify with the selection rules 
obeyed by the determinant interactions of fermions mediated by the classical instanton so- 
lutions [31, 32] of the non-Abelian gauge theory factors. More importantly, these conditions 
are independent of the Zn charge normalization, in contrast to the less physically moti- 
vated non-linear anomalies, Az? and Az$, which thus depend on the spectrum of massive 
decoupled modes. Among the linear anomalies, the gravitational and Abelian gauge U(l)y 
anomalies, A gra v and A±, are believed to be less robust than the non-Abelian ones, A3, Ai- 
Indeed, A\ is sensitive to the normalization of the hypercharge which remains a free param- 
eter as long as one is not concerned with the gauge group unification. Also, A gra v, Az 3 are 
sensitive to contributions from additional gauge singlet fields, which could either belong to 
the observable sector, such as the singlet S, or to the hidden sector, and hence coupled only 
through the gravitational interactions. While the non-linear anomalies A^2, Az'a are both 
sensitive to the ambiguity which arises due to the arbitrary charge normalization, Az* 
is also sensitive to presence of gauge singlets. 

We now discuss our results. These are displayed in Table I for the four different re- 
alizations of the three generalized parity symmetries in terms of the generator indices 
(m,n,p). Let us start first with the ordinary symmetries O. While it proves impossible 
to solve the complete set of equations, as already indicated, solutions do arise in very 
large numbers if one chooses to cancel the non-Abelian gauge anomalies .A3, A2 only. At 
this point, we mention that in the case of ordinary symmetries, and only for this case, 
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the cubic anomaly Aza automatically cancels once the equations for A3, A2 are satis- 
fied. The option of cancelling only the mixed gauge anomalies, ^.3, A2, Ai, (including 
the Az?> for the O symmetries is much more restrictive, but still yields solutions. The 
first realization of O symmetry occurs at the group order N = 9 with 3 GBP and GLP 
solutions and some 22 GMP solutions. We have quoted in Table I the values for the 
gravitational anomalies and Az? which remain generally uncanceled. For instance, the 
GBP generator (m,n,p) = (1,3,4) mod (9), has A gra v = —81 = mod (9), and hence 
a single uncanceled anomaly, Az* = —87 = —6 mod (9). At higher group orders, a re- 
stricted number of similar O solutions arise at orders N = 18 and N = 27. For the 
GMP case we have displayed in the table only a subset of the solutions. In fact, the 
O GMP solutions for the Z 9 group can be grouped into the three families of generators, 
(to, 3, 4), [to = 0,2,3,4,5,7]; (to, 6, 2), [to = 0, 1, 2, 4, 6, 7]; (to, 6, 8), [to = 0, 1, 4, 5, 6, 7]. 

Having a fixed N solution is not enough, since we must still solve the equations for the 
S field charges. A large number of solutions exist, in general, for the O symmetries under 
discussion. For instance, with the GBP solution (m,n,p) = (1,3,4) mod(9), the equations 
x+3y+4z = —3 mod(9), 3(x+3y+4z) = mod 9 admit about 80 different solutions of which 
we quote an illustrative sample: (x, y, z) — (0, 0, 6), (0,1,3), (0, 2, 0), (0, 3, 6), (0,4,3). The 
two other solutions, namely {m,n,p) = (5,6,2), (2,6,8) mod 9, have similar features. 

The case of Green-Schwarz anomalous symmetries O/GS is more constrained than the 
non-anomalous case discussed above. Unless one excludes the model dependent anoma- 
lies Az2, Az-i and A gra v, no solutions exist. However, by restricting again to the mixed 
gauge anomaly cancellation conditions only, which we express in terms of vanishing N 
linear combinations modulo N, As/k 3 — Aijki, kiAz/k^ — Ai, k\Aijk<i — A\, and set- 
ting the normalization parameters for the SM gauge group factors at the rational values, 
&3 = &2 = 1, k± — |, we find a single GBP solution appearing first for the group Z 7 with 
the generator {m,n,p) = (0,6,2). While this generator turns out to have a vanishing grav- 
itational anomaly, kgravAs/k^ — A gra v = —126 = mod 7, it still exhibits an uncanceled 
non-linear anomaly, Az? = —108 = —3 mod 7. Note that the next group order at which so- 
lutions appear is the integer multiple of the above with N = 14. Unfortunately, proceeding 
to the next stage of solving for the S field charges, we find it impossible to solve the relevant 
equations for the integers x, y, z. It is possible that this feature may be cured by adding an 
extra gauge singlet field. 
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The cyclic R parity discrete symmetries, R and R/GS, are more severely constrained 
than the ordinary ones. As seen on Table I, unique solutions are found at orders 12 and 
11, respectively, if one chooses to cancel the gauge anomalies only, while leaving A grav and 
the non-linear anomalies uncancelled. The R GMP solutions for the group Z 12 include the 
family (m, 6, 6), [m = 0, 1, 3, 4, 5, 7, 7, 9, 10], and the R/GS GMP solutions for the group Zj ii 
include the family (m, 6, 2), [m — 0, 1, 2, 4, 6, 7, 8, 9]. For both anomaly free and anomalous 
cases, we fail to find solutions for the equations for the S field charges x, y, z. 

In closing the discussion of results, we note that practically all the GMP solutions in 
Table I forbid the dimension 5 dangerous operators QQQL and U C D C U C E C . This is easily es- 
tablished by noting the corresponding selection rules, which require for O and R symmetries 
vanishing values for the total charges, g(QQQL) = — (n + p), g(U c D c U c E c ) = — (n — p) and 
A(QQQL) = -(2 + n+p), A(U C D C U C E C ) = -(2 + n-p), where A($ M ) = g(W)-(2-M). 

It is of interest to find out whether, by slightly modifying our search strategy, alternative 
options could exist. One might first consider solving the anomaly cancellation equations by 
setting the number of generations at the smaller values, N g = 1 or N g = 2. These solutions 
can be combined into quark and lepton generation dependent direct products, (Zn)% -± or 
(Zjv)tv 9=2 x (Zn)n 3= i. This option is not promising, however, since the solutions for N g < 2 
are even more scarce than for N g — 3. Thus, for GBP with N g = 2, the first solution for 
O symmetries occurs at N = 16, with (m,n,p) = (14, 12,6), for R symmetries at N = 12, 
with (m, n,p) = (8, 6, 6), and for R/GS symmetries at iV = 5 with (m, n,p) = (2, 4, 2). The 
case with N g — 1 does not have any solutions. 

Another option consists in enforcing the modulo N cancellation of A\ by adjusting the 
hypercharge normalization. As already noted, the freedom gained in relaxing A± anomaly 
cancellation constraint vastly increases the space of solutions. Changing the hypercharge 
normalization Y — > cY induces the modifications A\ — > c 2 A\ =^ k x — > k[ = c 2 k\. 
Specifically, given an O charge generator g with an uncanceled hypercharge anomaly, A\ ^ 
mod(iV), we can salvage the situation by transforming Y — > cY so that c 2 Ai = mod(A r ). 
For the anomalous GS symmetry, the same reasoning applies to the linear combination 
(?A\ — 2k[ = mod N. Of course, within a non-minimal grand unification or a string theory 
model in which some freedom is left for the hypercharge normalization, one must still consider 
how well the asymptotic prediction for the weak angle parameter, sin 2 9\y — ^2/(^2 + k[), 
fits in with the observed value. Considering, for concreteness, the illustrative case where 
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Ai = —pN + v, [p, v E Z + ] and changing Y — > cY, so that A[ = c 2 A\ = —pN, fixes the 
rescaling factor as c 2 = ^ = Sj^ +U — 1 + an d hence the modified asymptotic value of 
the weak angle as sm 2 6w ~ , where we have assumed, for simplicity, v « N. 

The case of an uncanceled A gra v anomaly may be treated by adding observable or hid- 
den sector chiral supermultiplet singlets, as already hinted above. Both kinds of singlets 
affect only A grav , Az$- Thus, given some generator with uncanceled gravitational and chiral 
anomalies, A gra v 7^ mod N, Azs 7^ mod N, one can attempt to rescue this solution 
by including an extra hidden sector singlet Si with R,A,L charges Xi,yi,Zi, and solving 
the equations, A gra v + Si = mod N, Az$ + Sf = mod iV with S = mx + ny + pz. 
Obvious modifications hold for the Green-Schwarz (GS) anomalous symmetries and for R 
symmetries. 

The option of extending the matter field content of the low energy theory by adding 
extra vector multiplets is also very efficient in relaxing the anomaly constraints. Indeed, 
we recall that the string theory models achieve consistency thanks to the presence of extra 
charged or singlet modes in the massless particle spectrum. To conclude, we note that 
interesting generalizations of our discussion would be to consider direct product of cyclic 
groups, Z N x Z M , or lepton flavor dependent cyclic groups. 



III. PHENOMENOLOGICAL IMPLICATIONS AND NEUTRINO MASSES 
A. Tree level neutrino masses 

After spontaneous breaking of the SU(2) L x U(1) Y gauge symmetry via the vacuum 
expectation values of the scalar components of H u , Hd, and S, the gauginos and Higgsinos 
mix with neutrinos. The resulting lepton number violating neutrino-gaugino-Higgsino mass 
matrix receives contributions from gauge interactions and the superpotential 

W v = \SH d H u + \ a L a H u S - (III.l) 

which arises from the last two terms of (II. 1) and the first term of (II. 2), respectively. The 
resulting mass terms of the neutrino-gaugino-Higgsino system can be written as 

C-mass = [\xH u H d + \v u H d S + \v d H u S - kxSS + H.C.] 

+ ^=-[v d H d - v u H u + H.c] - l -^-[v d H d - v u H u + H.c] 
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+ H ^a[xv a H u + v u v a S + v a H u S + H.c] 



(III.2) 



where A 3 is the third component of the SU{2)l gaugino X w , and A' is the U(l) gaugino. In 
(III. 2) we have used the following notation 



v u = < Hi >, v d =< H° d >, tan/3 = — , 



x = < S >, v a =< v a >, 



(III.3) 



while rest of the symbols have their usual meaning. Using (III. 2), we find the resulting 8x8 
neutralino- neutrino mass matrix in the field basis (—iX 1 , — iX 3 , H u , S, H d ^v e ^v^v T ) as 



M 



N 



Mi 





giv u 

V2 





-givd 

V2 


-givi 

V2 


-giV2 

V2 


-giv 3 

V2 





M 2 


-g2v u 

V2 





92V d 
V2 


92V1 
V2 


92 V 2 
V2 


92 V 3 
V2 


31 v u 


-giv u 





Y 


—Xx 


—X\x 


-X 2 x 


-X 3 x 


V2 










Y 


2kx 


Xv u 


-Xiv u 


-X 2 v u 


-~X 3 v u 


-givd 

V2 


92V d 
V2 


—Xx 


-Xv u 














~9l vi 
V2 


92V1 
V2 


—X\x 


-Xiv u 














-91V2 


92V2 
V2 


-X 2 x 


-~X 2 v u 














-giv 3 

I V2 


92 V3 
V2 


-X 3 x 


-X 3 v u 















where 

Y = -Xv d -^X a v a . 

a 

The mass matrix (III.4) can be written in the form 



M N = 



m 3x3 



where 



M x = 



Mi 








M 2 


giv u 


-g2v u 


V2 










-givd 


92 Vd 


L V2 


V2 



g\v u 

V2 





-gwd 


~92Vu 
V2 





g2Vd 





Y 


—Xx 


Y 


2kx 


-Xv u 


—Xx 


-Xv u 






(III.4) 



(III.5) 



(III.6) 



(III.7) 
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and 



m 



-givi 

V2 


-giV2 

V2 


-giv 3 

V2 


92 vi 
V2 


g2V2 

V2 


92 V '3 
V2 


— XlX 


— X2X 


— \ 3 X 


-Xiv u 
















(1IU 



The block form displayed in (III. 6) clearly demonstrates the "see-saw structure" of the mass 
matrix (III. 4). We note that the block matrix m characterizes the lepton number violation 
in the model. Furthermore, the NMSSM with lepton number violation is invariant under the 
S77(4) group acting on the set of superfields (Hd,Li), in the sense that the action of SU(4) 
transformations on (Ha, Li) leaves the superpotential form invariant up to corresponding 
transformations of Yukawa couplings. We can use this freedom to choose a basis which is 
characterised by vanishing sneutrinos vacuum expectation values, v a = 0. We shall choose 
such a basis in the following whenever it is convenient. 

The masses of the neutralinos and neutrinos can be obtained by the diagonalization of 
the mass matrix (III. 4) 

A/"*M N A/" _1 = diag(m x o,m Uj ) (III.9) 

where m x q,(i = 1, ....5) are the neutralino masses, and m u .,(j = 1,2,3) are the neutrino 
masses, respectively. The matrix (III. 4) cannot, in general, be diagonalized analytically. 
However, we are interested in the case where the tree level neutrino masses as determined 
from the mass matrix (III.4) are small. In this case we can find approximate analytical 
expression for the neutrino masses which are valid in the limit of small lepton number 
violating couplings. To do so, we define the matrix [38] 



If all the elements of this matrix are small, i. e. 



111.10) 



(n1.1i) 



then we can use it as an expansion parameter for finding an approximate solution for the 
mixing matrix Af. Calculating the matrix elements of ^ we find 

g 1 M 2 {\x) 



in 



V2det(M x o) 



[(2kx)x - 2Yv u ]Ai, 
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£i2 = 
£i3 = 
64 = 



g 2 M 1 (\x) 



V2det(A^ x o 
{g 2 M x + g\M 2 
2det(A4 x o) 



[(2kx)x - 2Yv u }A h 
[(2nx)v d x + (Xvl - Yv d )v u ]Ai, 



(gfMi + g 2 M 2 ) 
2det(.M x <>) 



{\v 2 u + Yv d )xAi, 



X lX (g 2 M 1 + g 2 M 2 ) 
65 " " A^ [1 + 2det(A< x o) (Av « " YVd) ] 
(g^Mi + g\M 2 ) 



4 2det(.M x o) 
where we have used the notation 



[(2Kx)v u xAi + (Xivi - Yvi)(Xvi - Yv d )\, 



(111.12) 



Aj = Xvi - XiV d . 



(111.13) 



We note that ^5 is not proportional to Aj. From Eqs. (III. 12) and (III. 13), we see that 
£ = in the MSSM limit where Aj = 0, v i = 0. The matrix A/"* which diagonalizes the 
neutralino-neutrino mass matrix M N can now be written as 



Af = 




V 



§«7 



(111.14) 



where we have retained only the leading order terms in £. The second matrix in (III. 14) 
block diagonalizes the neutralino-neutrino mass matrix M^r to the form diag (M. x o,m e ff), 
with 



m eff 



-m ■ Mj ■ m 1 



(M ig 2 + M 2 g' 2 )(Kx 2 -Yv u )x 
det(M x o) 



( A 2 e A e A u A e A\ 



A e A^ Al A M A 



\A e A T A^A T A 2 



;m.i5) 



where the tree level contribution to the light neutrino mass matrix admits the Feynman 
diagram representation of Figure 1(A). The eigenvalues of m e ff give the tree level neutrino 
masses. These eigenvalues are 

(M ig 2 + M 2 g' 2 )(Kx 2 -Yv u )x 



det(M x o) 



m 



m 



V2 



0, 



(111.16) 
(111.17) 
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(A) (B) (C) 

FIG. 1: The two point amplitudes for the process Vi — ► is? associated with contributions to the 
neutrino Majorana mass matrix at tree level (graph (A)) due to exchange of neutralino xu and at 
one- loop level (graphs (B) and (C)) due to intermediate propagation of neutralino xu sneutrino 
i)i and Higgs boson sector mass basis modes Si, Pi. The amplitude in (A) is initiated by tadpoles 
(VEVs) of H u , S, 9, and the amplitudes in (B) and (C) by double and single tadpoles of H u , Hd, S, 9. 
The cross on the neutralino propagator indicates a mass insertion term. 



where we have used m V3 > m V2 > m Vl . Thus, at the tree level only one neutrino is mas- 
sive. Its mass is proportional to the lepton number (and imparity) violating parameter 
X)jA?. A single lepton number violating coupling \ can lead to a non-zero neutrino mass. 
Furthermore, we note that the sub-matrices N and V v diagonalize M x o and m e ff. 



N*M x oN^ = diag(m x o), 
v J m effV„ = diag(0,0, ). 



(111.18) 
(111.19) 

Since only one of the neutrinos obtains mass, we can rotate away one of the three angles in 
the matrix V v . We can then write V v as a product of two matrices [39] 

( I \ ( cos9 13 " -oir,fl._ \ 

cos ^3 -sin ^3 ! 

sin #23 COS #23 J 



V 



where the mixing angles can be written in terms of Aj as 

tan# 13 = — T , 

tan 6*23 = -j^- 



sin 9 13 cos #13 j 



(111.20) 



(111.21) 
(111.22) 
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Finally, using the expression 



det{M x o) = (2Kx)(Xx)[(g 2 1 M 2 + g\M x )v u v d - M x M 2 {\x) 2 \ + 2M 1 M 2 {\v u ){\x)Y 



+ liglM.+glM^Xvl-vaYf, 



(111.23) 



for the determinant of A4 x o in (III. 16), we have for the tree level mass of the neutrino 

cos 2 (3 



m 



(111.24) 



where we have assumed that all the relevant masses (and the relevant vacuum expectation 
values) are at the electroweak (or supersymmetry breaking scale) scale denoted by in. For 
simplicity we have chosen the basis in which the sneutrino vacuum expectation values v i = 
to write the result (III. 24). We, thus, see that apart from the i?-parity violating parameter 
J2i{^i x ) 2 , the tree level neutrino mass is proportional to cos 2 (3. For large values of tan/3, 
this could lead to a suppression of m U3 , which could be important. 

It is now important to calculate the admixture of the singlet component (arising from 
the fermionic component of the Higgs singlet superfield S) in the three light neutrino states. 
From (III. 14) we can write the matrix JV* which diagonalizes the neutralino- neutrino mass 
matrix as 

-vji VJ(1 

The eigenvectors of the neutralino-neutrino mass matrix are then given by 



V 



(111.25) 



(111.26) 



where as indicated above we use the basis ipj = —i\ 3 , H u , S, Hd, v e , ^ v t)- The singlet 

component in the three neutrino states is then given by | A/64 1 2 , IA/74 1 2 and | A/s4 1 2 , respectively. 
For calculating these components we require the sub-matrix Vj£ of the matrix (III. 25). It 
is straightforward to calculate this sub-matrix, and the result is (A = (A e , A M , A r )) 

^ h N 

e 2 , (111.27) 

^ai|A| a 2 |A| a 3 |A| a 4 |A| e 3 + a 5 |A| y 



where 



ei = 



g e (A 2 + A 2 ) - A e (A^ + A r c T ) 



(111.28) 
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and 



= e A r + erA 
y/W + A?) 

63 = "v/(A e 2 + A^ + A?)' (m ' 30) 
s - ^ [1+ ( g lM 1 +^M 2 ) (A ^_ y ^ 



Ax 1 2det(.M x o) 

+ 2det(.M x o) [(A ^ - - YVd)l (IIL31) 

^ r (glM 1 + glM 2 ) 2 2 
C " = 2det(.M x o) ~ d 

+ 2det(.M x o) [{X ^ - YV ^ XV - - YVd)] ' (IIL32) 
_ _ \ 3 x. (^Mi + <? 2 M 2 V 2 v x2l 



Arc 2det(.M x o 

^^ lM - YM - Yv ' )] - (IIL33) 

(111.34) 



— -JSS)^- 2 ^ 1 ' (IIU5) 

a ^MS^) i(2Kx)x - 2Yv ' 1 (iil36) 

(g|Mi + £ 2 M 2 ) 2 
° 3 = 2 det(A4 o) K 2 * 3 ^ + ( At; " ~ Yv d)vu], (111.37) 

From (III. 25) and (III. 27) we obtain the important result 

IV 64 | 2 = |AT 74 | 2 = 0, (HI.40) 
|Ay 2 = a 2 |A| 2 . (111.41) 

Thus, at the tree level two light neutrinos do not have a singlet component, whereas the 
heaviest neutrino has a singlet component with a strength proportional to the square of the 
lepton number violating parameter |A|. 
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B. One-loop supersymmetry breaking contributions 

As shown above, at the tree level only one of the neutrinos obtains a mass through the 
lepton number violating Yukawa coupling A,, so that the tree level neutrino mass matrix 
can be written as 

ml = Vjm eff V v = diag(0, 0, m Ki ). (111.42) 

However, the neutrino mass matrix can receive contributions from loop effects. The super- 
symmetry breaking parameters are expected to play a crucial role through the one-loop cor- 
rections involving gauge interactions with exchange of sneutrinos [17, 20]. At one-loop level 
the needed suppression of neutrino masses can arise from cancellations between contribu- 
tions involving the Higgs sector, and from possible mass degeneracies among the sneutrinos. 
In the context of MSSM, this has been discussed in [24]. 

At one-loop level, finite Majorana neutrino masses can be generated through two classes 
of mechanisms involving either the gauge or superpotential interactions in combination with 
the soft supersymmetry breaking interactions. These mechanisms have been discussed in 
detail for the MSSM [29]. The loop amplitudes in the former class propagate matter parti- 
cles and contribute at orders AjAj, \\ijk and AjA' ijfc , and those in the latter class propagate 
sleptons and gauginos and contribute at orders A^.XiA^.Xj, A^.XiXj. These are associated 
with the mixing of sneutrinos and Higgs bosons, and are also responsible for the mass split- 
tings between sneutrinos and antisneutrinos [17]. The possibility that the combined tree and 
one-loop contributions in the MSSM could account for the observed flavor hierarchies in the 
masses and mixing angles of light neutrinos has been studied in several recent works [20-25] . 
The supersymmetry breaking interactions, the cancellations between contributions involv- 
ing the Higgs sector modes, and mass splittings among sneutrinos of different flavors, are 
expected to play a crucial role. 

The present section is aimed at studying the one-loop contributions to the neutrino mass 
matrix, and the extent to which these constitute a sensitive probe of the Higgs boson sec- 
tor of the NMSSM. The finite VEVs for the components of the scalar fields H d: H U: S, i>i 
can result in one- loop contributions for the two-point amplitude represented by the Majo- 
rana mass term L EFF = —\{m v )ijV C jV i + H. c. These contributions are displayed by the 
Feynman diagrams in Fig. 1(B) and (C), with double and single mass sneutrino-scalar mass 
mixing insertion terms. These are the analogs of MSSM [17] for the case of NMSSM. It is 
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important to carefully treat these see-saw like contributions by expressing the intermediate 
scalar and pseudoscalar propagators in the mass eigenbasis. We shall first obtain the scalar 
potential for the sneutrinos and Higgs bosons, minimize it with respect to the corresponding 
VEVs, va,v u ,x,Vi, extract the squared mass matrix whose off-diagonal blocks represent the 
sneutrino-Higgs mass mixing terms, and finally evaluate the one-loop contributions to the 
neutrino Majorana mass terms, (m u )ij. 



1. Coupling of Higgs boson and sneutrino sectors 



Using the standard procedure, we can write down the scalar potential of the NMSSM 
involving the relevant components of the complex scalar fields in terms of the F-terms, 
arising from the superpotential, the P-terms, arising from the gauge interactions, and the 
soft supersymmetry breaking terms as follows: 





= \H d H u S + XiLiH u S - 


■ -s 3 , 

3 ' 


V 


= V F + V D + Vsoft, 




v F 




gl+ / 2 i\v u \ 2 M 2 ) 2 


V so ft 


AB 


|2 i 2 1 |2 

v u \ +m s \x\ 



■[A- Xa \ a v a Vu + 



-x s + H. c.\. 



(111.43) 



In (III. 43), (pi stand for all the relevant scalar fields, and we have used the convenient four- 
vector notations A^ = (A, Aj), v A = (vd,Vi), with the summation convention over repeated 
indices undestood. The VEVs of different fields are temporarily extended to complex num- 
bers, v d =< H d >= v dl +iv d2 , v u =< H u >= v ul + iv u2 , x =< S >= x t +ix 2 , v { =< z>; >= 
Vn + ivi2, corresponding to the decomposition of Higgs boson and sneutrino fields into real 
scalar CP-even and imaginary pseudoscalar CP-odd components [3] 

vm + mi 



tj H dR + iH dI 
ti d — 1= , ri„ 



H uR + iH uI Sr + iSi _ 

S — 7= , Vi 



(111.44) 



V2 ' u V2 V2 V2 

This basis for the scalar fields (HdR, H u r, Sr), (H d i, H u i, Si) is related to the mass eigenstate 

basis (Si), (Pi), [I = 1,2,3] by the linear transformations 



(111.45) 



(H dR \ 












(H d i\ 








(PA 


H U R 


- 


V u 


= uj 


s 2 




H u i 


- V2% 


V u 




p 2 


V S R J 




\ X J 




\s 3 J 




\ Si J 




V x J 
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where U s , U p denote the unitary matrices which diagonalize the mass squared matrices in 
the interaction basis, (M 2 ) HdR>HuR>SR = M 2 sij and (M 2 ) HduHuuSl = M 2 tij , using the defini- 
tion U^ p M 2 p U s , p = (M 2 p ) diag . We then minimize the scalar potential with respect to the 
VEVs of various fields, and eliminate the soft supersymmetry breaking mass parameters 



s , mjj d M~.~.v* through the equations 

dV dV dV 1 , dV dV , . , 

— = 0, — = 0, -— = - i—) = 0. (111-46) 

dv u dx dv A 2 dv AR dv AI 

The mass squared matrices for the CP-even and CP-odd sector fields (H d , H U: S, h)R,i are 
then evaluated by applying the definitions 

M 2 .. = ^ , M 2 .. = £V . (111.47) 

Finally, we restrict our considerations to the physical vacuum solutions with vanishing imag- 
inary parts of the field VEVs, $>(v u ) = 0, $>(vd) = 0, ^s(x) = 0, $>(vi) = 0. For conve- 
nience, and without loss of generality, we shall also specialize to the choice of L A field basis 
characterized by vanishing complex sneutrino VEVs, Vi = 0. While feasible, the basis in- 
dependent analysis in the supersymmetry breaking case [19, 20] is significantly complicated 
by the need to account for several independent algebraic invariants. As a function of the di- 
mensionless parameters (A, X\, k), of the dimensional supersymmetry breaking parameters 
(A x , Ay, A K ), which include the gravitino mass parameter m 3 / 2 , and of the scalar field VEVs 
(vd, v u , x), the scalar and pseudoscalar mass squared matrices M 2 ^, M pi j, [i,j = d, u, S, Uk\ 
are given by the symmetric matrices 



Kdd = —\ 92 + 9l v z d + v u x(A\\ + k\x)], M 2 s4u = - 92 \ 9l VdV u + 2A 2 ^ U - A x \x - kX 2 x 2 , 
Vd £ £ 



M'i d s = -(A x \v u ) + 2X(Xv d - nv u )x, M 2 d0 . = —{A~ X X + kKx), 

^d 



Kuu = -[^r^< + v d x(A x \ + kXx)], M 2 m5 = -A x \v d + 2[- K \v d + (A 2 + X 2 )v u \x, 
M luui = 2XXiV d v u - x(A~ x .Xi + kXix), 

2 = A x Xv d v u + x ^_ AkK + 4k 2 x ^ M 2 = _ Axi \. Vu + 2X t (Xv d - kv u )x, 

X 

= M 2 ^ + &±£{& - v 2 u ) + XM + x% (111.48) 
for the CP-even scalars, and the symmetric matrices 

Kdd = —(A X X + k\x), M 2 p4u = x{A x X + kXx), 

Vd 
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Kds = Vu(A x \ - 2kXx), M 2 p4i>i = —{A~ Xi \ + K \ iX ), 

Kuu = —(A x \ + kXx), Ml uS = v d {A x \ - 2k\x), M 2 ^ = x{A x \ + k\x), 
M lss = 4K\v d v u + — — — + 3A k kx, M* SD . = v u (A- x \ - 2n\x), 

Ob 

K** = M h + 9l i^(vl - vl) + \M + x 2 ), (HI.49) 

for the CP-odd scalars. The orthogonal linear combinations of CP-odd scalar fields, G°(x) = 
cos f3Hdi{x) — sin f3H u i(x), A(x) = sin f3H,n(x) + cos (3H u j(x), identified with the decoupled 
Goldstone field which is absorbed as the longitudinal polarization mode of Z° and with 
the axionic symmetry pseudo-Goldstone boson mode A, respectively. The mass squared 
matrix in the field basis [A, Si, Ui\ is obtained by first applying the similarity transformation 
(G°, A, S/) T = Tl T (H d i, H u i, S 7 ) T , [K = diag(TZ p , 1)] with Hp denoting the SO (2) rotation 
of angle (3, and next by removing the decoupled Goldstone mode G°. The mass squared 
matrix in the transformed basis, (M 2 )g,a,5 7 = Ti r (M 2 )H dI ,H UI ,s I T^, = diag(1Zp, I4)] can 
then be written as 

M 2 pAA = sin 2 (3M 2 4d + cos 2 (3Ml uu + 2 sin (3 cos (3M 2 : ud = —1—^ X (A X \ + k\x), 
Kas = ^PKds + cos mlus = v(A x \ - 2k\x), 

Ml Ah = -^(A~ x ~\ t + K\ t x). (111.50) 

These results, with finite \x, are a generalization of the results of NMSSM [3] to the case 
when there is lepton number violation induced by trilinear couplings. These results reduce 
to the corresponding results of MSSM with lepton number violation [17] in the limit x — > 00 
with fixed \x = — /i, \x = —fa, and kx 3 . 

Since the off-diagonal entries in the sneutrino-Higgs mass matrix can be safely assumed 
to be small in comparison to the diagonal entries, one may evaluate the contributions to the 
sneutrino mass splittings by making use of second order matrix perturbation theory. The 
same approximation is also used in evaluating the modified sneutrino propagators in the 
mass insertion approximation. Specifically, the Higgs boson sector propagator, modified by 
tbe two mass mixing terms v, L x {Si © P/) x i/j in the Feynman diagram of Fig. 1(B), can be 
represented by the weighted propagator 

*W«) - E 3^r- E en-si) 

J=l,2,3 q M Sj J=l,2 q Mp J 
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where 



Pij 



E M?*k(vT)kj)( E <^(^ T w), (HI.52) 

- J - c 7 K k=d,u,S 7 

(111.53) 



k=d,u,S 



k=d,u,S k=d,u,S 



Evaluating the transition amplitude for the double mass insertion one-loop Feynman graph 
of Fig. 1(B) with the above formula for the weighted Higgs boson propagator, one obtains 
the contribution to the light neutrino mass matrix 



= jE M x ; (^2-tan^iV n ) 2 [ J2 sj 3 h(m^m,M^M Sj ) 



J=l,2,3 



E Pijh K , m p . , il% , M Pj )\ , 



(111.54) 
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[/ 4 K,m s .,M^,M Xj ) 
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d 4 q 



Cim^mo^M^Mxj) 



z(2tt) 4 (g 2 - m\)(q 2 ~ m%)(q 2 ~ M? )(q* - M 2 J J ' 



(111.55) 

where tan0 w = g\/g2, and we have used the matrix iV to denote the unitary transformation 
linking the interaction and mass eigenstates of massive neutralinos, (xm)mass = (N*) m ixi- 
The momentum integral I4 admits the analytic representation [40] 

I 4 (m 1 ,m2,m 3 ,m 4 ) = — ^ — f [h(m 1 ,m 2 , m 3 ) - / 3 (mi, m 2 , m 4 )], 

77^3 — 777.4 



73(7771,7772,777,3) 
7 2 (r77i,777 2 ) 



2-[ J 2(^l,^2) - I 3 (^l,^3)], 

7772 - m 3 



1 



7/1 



-loe — 

(47f) 2 7772 ~~ m l m 2 
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The single mass insertion one- loop amplitude, displayed in Feynman graph of Fig. 1(C), 
yields the following contribution to the light neutrino mass matrix 



4 , My. 

Z,m X; 



iV w (A^ 2 - tan0 w JV ml ) + N m4 (N l2 - tan 6 w N n ] 



x 



[E(^W^,-E(^W3y 



j=i 



j=i 



-iV i3 (iV m2 - ixm6 w N„ a )[E,(Vj')ujQl i j ~ Y,K)ujQ p ^ 

j=i j=i 



+ (i <-» j), 
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where various quantities in the above equation are defined as 

h{m 0v M^M Xj ) 



An examination of the off-diagonal matrix elements of the sneutrino-Higgs squared mass 
matrix given by equations (III. 48) and (III. 49) shows that the above one-loop contributions 
to the neutrino mass matrix consist of sums of two separate matrices involving the three- 
vectors, \ and Ay\ in the space of fields Lj. Combining these with the tree contribution 
discussed in subsection III A, one can now write the following representation of the effective 
light neutrino mass matrix as a sum of three contributions 

(m„) y = (X^ + X^ + X^KXj + Y^KAx^ 

+(Z l B + Ztt&A^Xj + A^XiXj), (111.62) 

where the lower suffix labels A, and B, C in the coefficients X, Y, Z refer to the tree and 
one-loop contributions coming from the Feynman diagrams (A), and (B), (C) in Fig. 1 and 
we have appended the upper suffix labels t, I to emphasize the distinction between tree 
and one-loop contributions. We note the absence of the coefficient Yq, and the relation 
X l B ,X l c << X\ expected from the loop suppression factor, which allows us to ignore the 
coefficients Xb and X l c . The single mass insertion contributions Z l B , Z l c , and the double 
mass insertion term Y B have the ability, either separately or in combination, to produce 
a second non- vanishing mass eigenvalue, provided only that the three- vector A^.\ is not 
aligned with Aj. We recall that the three- vector proportionality, A^\ oc Aj, would hold if 
supersymmetry breaking were flavor universal. Moreover, as was first observed by Chun et 
al. [23] in the context of MSSM, application of matrix perturbation theory to the additively 
separable neutrino mass matrix (m u )ij = x/ii/ij + ybfij + z(^bj + fijbi) indicates that the 
two finite eigenvalues present in the limit y, z << x are given by xfij + lybiHi + 0(y 2 , z 2 ) 
and yb 2 + 0(y 2 ,z 2 ). Hence, assuming in the above Majorana neutrino mass matrix (m u )ij 
that the coefficients Y B , Z l BC are of subleading order relative to X l A , one concludes that 



E Ml 0jk {Uj) k jh{m^M^M Sj ), 

k=d,u,S 

E Ml- jk (Uj) kJ I 3 (rn 0i ,M^,M 2 Pj ), 
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the second non- vanishing eigenvalue is of first order in Y B but of second order in Z B , Z]j, 
namely, m V2 ~ 0(Y B )+0(Z l B , Zj§). Thus, as far as the second mass eigenvalue is concerned, 
this implies that the single mass insertion amplitude (C) is subdominant, so that we can 
restrict consideration to the double mass insertion contribution (B) only. 

A rough estimate of various contributions can now be obtained by isolating the stronger 
dependence on tan/3, while assuming that all the mass parameters take values of same 
order of magnitude as the supersymmetry breaking mass scale rho- This yields the approxi- 
mate formulas for the coefficients representing the tree and one-loop contributions, X A and 
Y B , Z l B , Z l c , respectively: 

t x 2 cos 2 f3 , x 2 

m m cos z p 

Z l B * z KX \n e L e H , Z l c ~ J^—e L e B , (111.63) 
m cos 2 p m cosp 

where we have included the suppression effect from the loop in the factor ~ l/(4ir) 2 ~ 
1CT 2 and that from the Higgs sector in the factors e B and e' H . Note that we have omit- 
ted the one- loop contributions to the component XiXj, which are associated with the sup- 
pressed coefficients, X l B ~ Z l B , X l c ~ Z l c . The Higgs sector decoupling effect arises from 
the cancellation between the contributions from CP-even and CP-odd scalars, and is most 
effective in the case where the lightest scalar contribution is well separated from the other 
modes, and the mass spectrum is ordered as, m Sl ~ m z « m S2 ,m S3 ,mp 1 ,mp 2 . The 
dominant contribution to the second finite neutrino mass eigenvalue is, then, of order 
m V2 ~< A\ 2 \ 2 > cos 2 PeLtM/mo- Moreover, a third finite mass eigenvalue may be generated 
from the one-loop amplitudes under study by taking into account the flavor non-degeneracy 
in the sneutrino mass spectrum. The presence of a small relative mass splitting for the sneu- 
trinos, say, ui, has the ability to produce a third non- zero neutrino mass eigenvalue. The 
relationship between the ratio of non-zero neutrino masses and the sneutrino mass splitting 
is given by [24] 

e D = ^ ^ lO^A 2 , 

m U2 

A e = K-<l . (IH.64) 
2m^ 
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2. Higgs boson decoupling 



Our next task is to estimate semi quantitatively the Higgs sector suppression factor e#. 
Following MSSM [24], we consider the definition 

I E i s\j h ijn Vi , m v . , My_ , M Sl ) - E j pj h {m* , m v . , , M Pj ) | 



en 



(111.65) 



Ej Wijhim^.m^, M Si ) \ + Ej \pijh(m Vi , m Vj , M%, M Pj )\ ' 

One can easily verify that e# vanishes in the limit of mass degenerate scalars and pseu- 
doscalars. The Higgs boson mass eigenvalues m 2 Sl (I = 1,2,3) and m 2 Pj , (J = 1,2) and 
mixing matrices U s and U p , where the latter is expressed in the basis (A, S) as U p = TZ 1 
in terms of the SO (2) rotation matrix of angle 7, are determined once one substitutes the 
values of the free parameters A, k, Ax, A k , tan/3 = v u /vd, while using the observed value 
v = (v'j + v 2 )^ = 174 GeV. Ensuring a vacuum solution with electroweak symmetry breaking 
at the appropriate scale, and without tachyonic scalar modes, is known to impose strong 
constraints on the NMSSM [3]. However, a systematic exploration of entire parameter space 
consistent with all the physical constraints is beyond the scope of the present work. For a 
semiquantitative estimate, which is adequate for the purpose of illustrating the typical order 
of magnitude values assumed by e#, we only explore a small region of parameter space. For 
this purpose, we shall consider a modest numerical study confined to large and small values 
of x, respectively, with the coupling constants held fixed, where one expects to find the 
largest departures from the MSSM. The mixing matrices in these two regimes are described 
by the approximate formulas 



x » vi, v 2 : U s ~ diag(Tl(3, 1) , 7 ~ vr/2, 

/ — sin a + C cos a cos a + C sin a \ 



x « Vi, v 2 : U s 



V 



cos a + C sin a sin a — C cos a 
I) 



(111.66) 
, 7 ~ 0, (111.67) 



where C = 2\A x xcos(2a) sin2( / 9 - a)/[m|sin(4/3)]. 

We shall study the dependence of the function e# on various parameters by means of 
two different prescriptions. In the first, we set the dimensionless coupling constants at 
the renormalization group infrared fixed point values, A = 0.87, k = 0.63, with the soft 
supersymmetry breaking parameter A K having the fixed value A K = 200 GeV, and vary 
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the parameter A\. This results in the variation of the mass m c of the charged Higgs 
boson, C + = cos f3H+ + sin (3H^*, which is given by the tree level formula m c = m^ — 
\ 2 {v1 + v 2 a ) + 2\(A X + kx)x/ sin(2/3). The variation of e H with m c is examined at discrete 
values of x and tan (3, while the sensitivity of these results with respect to the other fixed 
parameters is examined by considering small variations around the above reference values of 
the parameters. In the second prescription, we examine the dependence of e# on x for the 
choice of fixed parameter values A = 0.5, k = 0.5, A K = 100 GeV, A x = 100 GeV. In both 
prescriptions, we assign definite mass values to the lowest lying neutralino and the pair of 
lowest lying sneutrinos, namely, m^ = 300 GeV and m^ = 100 GeV, m„ 2 = 200 GeV. while 
noting that the loop momentum integral I4 depends very weakly on the input masses. 

The plots of the ratio e# as a function of mc and x is displayed in Fig. 2 in the frames 
(a), (6) and (c) for the above two prescriptions. For the first prescription using variable 
mc, the plots are restricted to the physically acceptable values of the parameter A x in 
which no tachyonic scalars or pseudoscalars are present in the neutral Higgs boson sector. 
In the small x regime with r = 0.1, the lowest lying Higgs boson mass lies in the interval 
mh ~ 50 — 20 GeV for mc ~ 20 — 100 GeV, which is excluded by the experimental 
limits. In the intermediate x regime with r = 1 and r = 10, it is pushed up to the interval 
m h ~ 70 - 120 GeV and m h ~ 130 - 140 GeV, respectively for m c ~ 100 - 300 GeV and 
m c ~ 300 — 2000 GeV. The plots in the frames (a), (b) show that the variation of e H with 
mc is slow except when one approaches the boundaries where tachyons appear. The typical 
size of the ratio is e# = O(10 _1 ), irrespective of the small or large values of x, but decreases 
by a factor 2 — 3 with increasing tan/5. The discontinuous behavior of the curves for en 
is explained by the fact that this is the absolute value of the difference of two amplitudes. 
The plot in frame (c) shows that en has a strong variation with increasing x in the interval 
x > v with a typical size O(10 _1 ), decreasing by a factor 10 with increasing tan/5. 

We have also examined how the ratio en varies with small variations about the reference 
values of the couplings for the first prescription. Changing A K has a mild influence on the 
Higgs boson mass spectrum and hence on tu- Indeed, increasing A K by a factor 2 — 3 does 
not affect the prediction for en significantly. Decreasing A by a factor 2 reduces en mildly 
at small x and more strongly, by factors of 2 — 5, at large values of the parameter x. A 
similar but weaker decrease applies when we reduce k by a factor 2. As shown by Fig. 2, 
en decreases rapidly with increasing tan /?, but undergoes very small changes when we allow 



29 



for large variations of m^, m^. 

Thus, the main conclusion of our analysis is that the suppression factor e# arising from 
the Higgs sector is typically of order 1CT 1 — 1CT 2 . This is larger than the value obtained 
for the corresponding factor in the MSSM [24], e# — 10~ 2 — 10~ 3 , at the values of param- 
eters consistent with physical constraints. It is, however, possible that there are regions of 
parameter space where e# is significantly smaller in the NMSSM. 



IV. FLAVOR SYMMETRIES 

There are too many parameters in supersymmetric models, including the nonminimal 
super symmetric model, to make any specific predictions for the neutrino spectrum. It is even 
difficult to identify important contributions to the neutrino masses. Here we shall study a 
specific framework, that of an Abelian flavor (horizontal) symmetry [41, 42], where specific 
predictions can be made. Flavor symmetries are usually invoked to explain the pattern of 
fermion masses. However, any theory of fermion masses must also explain why the violations 
of i?-parity (or lepton and baryon number) are small. This applies particularly to NMSSM 
with lepton number violation coming from a trilinear type superpotential coupling that we 
are considering here as the origin of neutrino masses. 

We start by recalling the salient fatures of the Abelian flavor symmetry framework. The 
basic idea is to use an Abelian horizontal symmetry U(l)p to forbid most of the Yukawa 
couplings except perhaps the third generation couplings. The hierarchies of fermion masses 
and mixing are then generated through higher dimensional operators involving one or more 
electroweak singlet scalar fields. These fields acquire vacuum expectation values at some 
high scale and give rise to the usual Yukawa couplings. More specifically, if 6 is some such 
field which has charge —1 under U(1)f, then X-charge allows the non-renormalizable term 
in the superpotential 

VWr(-^J , (iv.i) 

where $j is a matter superfield of flavor i, and H is a Higgs superfield with appropriate 
transformation properties under the gauge group. The coupling \j is of order unity, and 
M is some large mass scale. The positive rational numbers are nothing but the sum of 
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FIG. 2: The Higgs sector decoupling ratio en in the double mass insertion approximation at one- 
loop order is plotted as a function of the charged Higgs boson mass mc for the three regimes of 
the VEV ratio parameter, r = - = 0.1 and 1, 10 in the frames (a) and (b), respectively, and as 
a function of the S field VEV, x, in the frame (c). The results in frames (a), (b) are obtained 
with the fixed values of parameters A = 0.87, k = 0.63, A K = 200 GeV, and with a variable 
parameter A\, which determines the charged Higgs boson mass mc- The plots in frame (a) are 
for tan/3 = 1.5, 4, 10, with r = 0.1, and those in frame (b) for r = 1, tan/3 = 1.5, 4, and for 
r = 10, tan/3 = 1.5, as indicated in the legends. The curves in frame (c) are for tan/3 = 1.5, 4, 16 
at the fixed values of parameters, A = 0.5, k = 0.5, A K = 100 GeV, A\ = 100 GeV, as indicated 
in the legend. The variables in the 4-point amplitude I4 are chosen by setting the masses of the 
lightest pair of sneutrinos at mo 1 = 100 GeV, m^ 2 = 200 GeV, and that of the lightest neutralino 
at = 300 GeV. 
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n ij ~ V* T <fj 

When gets a vacuum expectation value, an effective Yukawa coupling 
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= Xij9 n c ij , (IV.3) 

is generated. If 9c is a small number, and if the U(1)f charges are sufficiently diverse, one 
can implement various hierarchies of fermion masses and mixing. This can then be viewed 
as an effective low energy theory that originates from the supersymmetric version of the 
Froggatt-Nielsen mechanism at higher energies. From above we then have the following 
consequences: 

(i) Terms in the superpotential that carry charge n > are suppressed by 0(9 C ), whereas 
those which have n < are forbidden by the holomorphy of the superpotential. 

(ii) Soft supersymmetry breaking terms that carry a charge n under U{l)x are suppressed 
by 0{9^). 

Applying the above scheme to the neutrino mass matrix, we see that the additive sep- 
arable structure of the combined tree and loop level contributions give us the ability to 
account for moderate flavor hierarchies. Let us first recall that the fit to the neutrino 
oscillation experimental data, assuming a mass spectrum of normal kind with mild hierar- 
chies, favors the following approximate solution for the three masses and mixing angles [26] : 
m V3 ~ 10 _1 eV, m U2 ~ 10~ 2 eV, m Ul ~ 10~ 3 eV and sin 2 9 2 3 ~ |, sin 2 9 12 ~ §, sin 2 9 13 < 
1 .4 1CT 2 . Of course, the contributions that we have discussed so far are controlled by 0(100) 
GeV weak interaction scale, which lies considerably higher than the observed neutrino mass 
scales. Having identified the supposedly dominant contributions, it is now necessary to find 
a plausible suppression mechanism which accounts for the wide O(10 12 ) hierarchy in mass 
scales. As in the familiar Froggatt-Nielsen approach [41, 42], we can adjust the overall size 
of the contributions without an excessive fine tuning of the free parameters by postulating 
that the superpotential and supersymmetry breaking couplings of the NMSSM arise from 
non-renormalizable operators with effective couplings weighted by powers of the small pa- 
rameter 9 C =< > /M, which we shall identify here with the Cabibbo angle parameter, 
Q c — 0.2. With h(Li), h(Hd jU ),--- denoting the Abelian horizontal group U(l)p charges 
assigned to the various superfields, one finds, \ = 0^( L ^ +h ^ +h( - H ^ < \. > an( j similarly 
for the associated supersymmetry breaking parameters, A^.Xi, with the expectation that 
< A 2 >= 0(1), < A| Af >= 0(1). Using the results of subsections III A and III B, one can 
write the predicted finite neutrino mass eigenvalues as 



x2 c °s 2 P ^ r 2 . a 2h(L 3 ) rn U3 cos 4 f3 < A 2 > 2 [h(L 3 )-h(L 2 )] 
< a > Vq , ~ Tjry— &c - 
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^^e D ef^- h ^ )] . (IV.4) 

In order to obtain m Ki ~ 10" 1 , we must have O 2 ^ 1 '^ > 1CT 12 , and hence h(L 3 ) < 9. Similarly, 
in order to obtain m U2 ~ 1CT 2 yields h(L 3 ) — h(L 2 ) < 2. Furthermore, m vi ~ 1CT 3 can be 
achieved by lifting the mass degeneracy between sneutrinos with h(L 2 ) ~ h(Li). Recalling 
the predictions for the lepton flavor mixing angles, sin^- ~ I ^ L% ^ h ^ Lj ^ it follows that, as 
in the case of MSSM case [24], the selection of horizontal symmetries involving nearly equal 
horizontal charges h(Lj) introduces a fine tuning problem in order to account for the small 
observed mixing angle #13. 



V. SUMMARY AND CONCLUSIONS 



We have studied the nonminimal supersymmetric standard model with lepton number 
violation in detail. This model has a unique trilinear lepton number violating term in 
its superpotential, and a corresponding soft susy breaking scalar trilinear coupling. We 
have attempted to justify these terms on the basis of a gauged discrete symmetry. We 
have shown that these terms give a viable description of the light neutrino Majorana mass 
matrix provided one stabilizes the large mass hierarchy with respect to the weak gauge 
interactions scale by invoking horizontal flavor symmetries. A satisfactory feature of this 
extended version of the NMSSM is that the suppressed interactions are all associated with 
effectively renormalizable and dimensionless Yukawa couplings. This mechanism represents 
an economic alternative option to the familiar see-saw mechanism of generating the light 
neutrino mass matrix. Although qualitatively similar to the bilinear lepton number violation 
that occurs in MSSM, it distinctly differs from it on important quantitative grounds. We 
find that only one of the three neutrinos obtains mass at the tree level and that this has 
a finite component of the massive singlet fermion. We have also calculated the one-loop 
radiative corrections to the neutrino mass matrix generated by the coupling of sneutrino 
and Higgs boson sectors. These can contribute finite masses to the other two neutrinos in a 
manner favoring mild hierarchies of normal kind for the neutrino mass spectrum along with 
large lepton flavor mixing angles. One can reproduce a single small mixing angle, as needed 
for agreement with the current experimental data, at the cost of a small fine tuning. 

In an effort to put on a firmer theoretical basis the different versions of the NMSSM 
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with renormalizable B or L number violation, we have also studied the four main gauged 
Zn cyclic group (ordinary and R, free and GS anomalous) realizations of the generalized 
baryon, lepton and matter parities. The constraints from the anomaly cancellation condi- 
tions are so strong that no solutions exist if one restricts to the strictly minimal matter 
field content. However, making the reasonable choice of retaining only the least model de- 
pendent conditions, associated with the mixed gauge anomalies, and of admitting at least 
one extra gauge singlet chiral superfield (in addition to the standard one which couples to 
the Higgs bosons), we find interesting restricted classes of symmetry solutions at low cyclic 
group orders N. 
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APPENDIX A: CLASSIFICATION OF CYCLIC DISCRETE GAUGE SYMME- 
TRIES 

The interest in generalized parities for the MSSM was historically motivated by the need 
to suppress the dimension-5 baryon and lepton number violating supersymmetric opera- 
tors [35, 36]. Our purpose in the present appendix is rather to classify the discrete cyclic 
group symmetries which protect the structure of renormalizable versions of the NMSSM 
superpotential with baryon or lepton number violation. The issue of adding gauge singlet 
chiral supermultiplets was considered by Lola and Ross [37], although their work was fo- 
cused on applications involving non-renormalizable couplings. Our present treatment of this 
problem also slightly deviates in certain technical details from that followed in this earlier 
work. 

We wish to prove the existence of cyclic symmetries of general order N which leave 
invariant the interaction superpotential of the NMSSM with B or L number violation with 
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trilinear Yukawa interactions. Based on the approach initiated by Ibanez and Ross [35, 36], 
one distinguishes three cases of discrete symmetries designated as generalized baryon (GBP), 
lepton (GLP) and matter parities (GMP), respectively. For each case, there are four dif- 
ferent realizations depending on whether the symmetry is ordinary or R like, and whether 
it is Green-Schwarz (GS) anomaly free or anomalous. We discuss first the general classi- 
fication of the different discrete symmetries and next the consistency conditions imposed 
by the cancellation of anomalies. Our considerations will be restricted to the flavor blind 
symmetries. 

1. Ordinary symmetries 

Let us start with the ordinary anomaly free symmetries. Recall first that the quark and 
lepton generation independent Abelian charges conserved by the renormalizable R parity 
conserving (RPC) superpotential couplings of the MSSM form a vector space generated by 
three continuous U(l) symmetries. A convenient basis for the three independent charges is 
given by R, A, L, where R = T 3R and A — Y A identify with the Cartan generators of the 
right symmetry group SU(2) R and the SU(2) group embedded in SU(6) x SU(2) C E 6 , 
and — L identifies with the usual lepton number. The charges R, A, L assigned to quarks, 
leptons and Higgs boson superfields are displayed in the following table, along with those 
assigned to the singlet superfield Si, which are denoted by x, y, z. (The charge P will appear 
in the next subsection in the discussion of R symmetries.) 
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Mode 


Q 


u c 


D c 


L 


E c 
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QY 


1 


-4 
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-3 
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-3 


3 





R 





-1 


1 
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-1 


1 




A 

1 1 


n 


n 






n 


i 

± 


n 


y 


L 











-1 


1 








z 


P 


-1 


-1 


-1 


-1 


-1 


1 


1 


-1 



The Z§, Zft, Z^ group elements are constructed in the same way as for the continuous 
groups, U(1)r : a,l, by writing R = e taRR , A = e taAA , L = e iaLL , while restricting the 
complex phase angles to the fixed values, an = 2nm/N, aA = 2nn/N, ai = 2np/N, with 
integer charges m,n,p defined modulo N. Note that the U(\) P q symmetry generated by 
gpQ = R 2 A is a chiral, color group anomalous symmetry which conserves all renormalizable 
(RPC and RPV) trilinear couplings of the MSSM. The pseudoscalar Higgs boson A is the 
pseudo-Goldstone boson of the U(1)pq symmetry present in the limit fi — > where the 
explicit symmetry breaking bilinear coupling \iH u Hd is absent. 

The multiplicative Z^ symmetries of the renormalizable superpotential for the quarks, 
leptons and Higgs bosons may be parameterized in terms of the generators g = R m A n L p = 
gpQR m ~ 2n L p , [m,n,p integers]. The symmetry solutions preserving B and L, or B 
alone or L alone are designated as generalized matter, baryon and lepton parities (GMP, 
GBP, GLP), respectively. Thus, aside from the regular interactions with Higgs bosons, 
QU C H U , QD c H d , LE c H d , the GBP are required to forbid the interactions U C D C D C but to 
allow the interactions LH U , LLE C , LQD C , the GLP acts in a manner forbidding the lepton 
number violating interactions, but allowing baryon number violating interactions U C D C D C , 
while the GMP must forbid all the matter interactions. Noting the charges for the pure 
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matter couplings, g(LLE°) = g(LQD c ) = m — 2n — p, g(U c D c D c ) = m — 2n, one finds that 
the discrete symmetry generators preserving the MSSM trilinear superpotential in the three 
relevant cases are given by [35, 36] 

GBP : m - 2n - p = 0; m - 2n ^ =S> g GBP = g n PQ {RLf , [p + 0] 

GLP : m - 2n = 0; m - 2n -p ^ 0, ^ =}► (? GLP = ^ Q L P , [p ^ 0] 

GMP : m - 2n - p ^ 0, m - 2n ^ 0, =>• gcMP = g PQ R m ~ 2n L p , [m-2n^0, p ^ 0]. 

(A.l) 

A similar analysis applies in the NMSSM, with the generators for GBP, GLP and GMP 
required to forbid the matter couplings violating baryon number only (U C D C D C ), lepton 
number only (LiH u S, LLE C , LQD C ) and both combined, respectively The selection rules 
for the allowed and forbidden S field dependent trilinear couplings are given by, g(H d H u S) = 
n + S = 0, g(LH u S) = m — n — p + S = 0, g(S 3 ) = 3S = 0, [5 = mx + ny + pz] and 
g(S') =5^0, ^(S* 2 ) = 25 ^ 0. Except for the different conditions on the integers 
(m,n,p) G Z N , the cyclic symmetry generators have the same functional form as in the 
MSSM, 

GBP : m-2n-p = 0, m-n-p + S = 0; m-2n^0 =>- go bp = g PQ (RL) p , 

GLP : m — 2n = 0; m — 2n — p 7^ 0, n — p + 5 ^ , m — n — p + iS^0 =>- #glp = 9 P qL p , 

GMP : m-2n-p^0, m-2n^0, m-n-p + S ^0 =^ g GMP = g^ Q R m - 2n L p , (A.2) 

where the three symmetry cases must satisfy p ^ along with the conditions: n + S = 
0, 3S = 0;n ^ 0,S ^ 0,25 ^ 0. 

Having classified the cyclic groups, we now wish to implement the condition that these 
belong to gauged symmetries. This means that the total contributions to the quantum 
anomalies from massless fermions of the low energy theory must either vanish or be com- 
pensated by those of the massive fermions of the high energy theory which decouple by 
acquiring large Dirac or Majorana masses. The coefficients of mixed gauge and gravita- 
tional anomaly operators F a F a , RR, FgFy and of chiral anomaly operator, Fjj, acquire the 
following contributions from massless fermions, 

A{Z N x G 2 a ) = Yl^i)g(^i), A(Z N x grav 2 ) = 

A(Z%) = E^ 3 (4), A(Z 2 N x U(1) Y ) = ££ 2 (^ (A.3) 

i i 
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where /i a (V , «) denotes the Dynkin index of the fermion ipi representation with respect to 
gauge group factors G a of the SM gauge group SU(3)c x SU(2) L x U(l)y and the label 
'grav' stands for the gravitational field source. The Abelian gauge anomaly Ai will be 
evaluated by setting conventionally the hypercharge normalization such that, Y(e c ) = 1, 
which implies that the trace of Y 2 over a single quark and lepton generation amounts to 
Trace(Y 2 ) = 2k\ = 10/3. The anomaly cancellation conditions for the NMSSM are given 
by the formulae 



A 3 = A{SU{3) 2 x Z N ) = -nN g = rN, 

A 2 = A{SU{2) 2 x Z N ) = -N g {n + p) + N 2h n = rN, 

A, = A(U(1) 2 Y x Z N ) = N g (-^ + |) + N 2h ^ = rN, 

A Z 2 = A(U(1) Y x Z 2 N ) = -2N g (2mn + p(n - m)) - N 2h n(n - 2m) = rN, 

N 

Ag rav = A(grav 2 x Z N ) = -N g (5n - m + p) + 2N 2h n + S = rN + i]s— , 
A z z =A(Z%) = N g [-3(m 3 + (n - m) 3 ) - 2{n + pf + (m + pf\ 



Here, N g denotes the number of quark and lepton generations, N 2 n the number of H u 
Higgs boson supermultiplet pairs, and the symbols r, s G Z denote arbitrary integers (taking 
independent values for the different anomalies) so that the different equations are understood 
to be satisfied modulo N. The additional vanishing conditions associated with the parameter 
i] = 0, 1 for iV odd and even, respectively, are introduced to account in the even N case for 
the presence of massive Majorana fermions in real representations of the gauge group factor 
G a . The gauge singlet charges enter only through the linear combination, S = mx + ny + pz, 
which is set to S = —n. 

It is straightforward to generalize the above results to the case involving several gauge 
singlet chiral supermultiplets, One just needs to assign Si the R, A, L charges Xi,yi,Zi, 
and to replace in the anomaly coefficients, S — > J2i<Si, <S 3 ~^ J2i<Sf. These additional 
contributions set conditions on the charges Xi,yi,Zi expressing the net cancellation of the 
anomaly coefficients A grav , Azs- Each allowed coupling must also be accompanied by an 
additional constraint equation expressing the associated selection rule. 



+N 2h [(n - mf + m 3 ] + S 3 




(A.4) 
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The following two-stage procedure may be used in solving the anomaly cancellation condi- 
tions for each fixed N generator. One first scans through the non- vanishing integers m, n, p 
to select those satisfying the above set of equations and next scans the non-vanishing inte- 
gers x,y,z G Z]\f which solve the equations S = mx + ny + pz = —n 7^ 0, IS ^ 0. The 
search is most easily implemented with the help of a numerical computer program. 

We study next the Green-Schwarz anomalous discrete symmetries. This case differs 
from the anomaly free one in that the anomaly coefficients in the effective action is now 
allowed to take finite values, provided only that these are canceled by the additive contri- 
butions to the anomalies of universal form associated with the gauge and gravitational cou- 
plings of the model independent axion-dilaton chiral sup er mult iplet. The modified anomaly 
cancellation conditions are then expressed in terms of the shifted anomaly coefficients, 
A a —2k a SGs — 0, Agrav—2kg rav 5Gs = 0, [a = 3, 2, 1] where k a are rational parameters (integer 
quantized for non-Abelian group factors G a ), k grav = 12, and Sqs denotes a universal model- 
dependent parameter reflecting the underlying high energy theory. These conditions can 
also be represented by the proportionality relations, Az/k^ = Aijki = Ai/ki = A grav lYl. 
The parameters k a in the minimal gauged unified theories are set at the numerical values, 
h = h = 1, h = 5/3. 

2. R symmetries 

We shall continue using the abbreviations GBP, GLP, GMP for the generalized baryon, 
lepton and matter R parity discrete symmetries. A convenient representation of the Z N 
group generators can be constructed by introducing the fermionic generator P defined by 
its action on the superspace differential, P ■ d6 = e~ 2m l N d9 and by the charge assignments 
of the gauginos, P(g) = P(W) = P{B) = 1, and of the matter and Higgs boson superfields, 
as displayed in the table placed at the beginning of Section A. With the understanding 
that the charge assignments displayed in the table for P and for R, A, L apply to the 
fermion field component of the chiral superfields, the Z N group generators preserving the 
MSSM matter-Higgs boson trilinear superpotential are simply given by, g = PR m A n L p The 
g charges of fermion and scalar field components of chiral superfields ip, <fi are then related 
in the usual way, g(ip) = g{<t>) — 1, so that a superpotential term W is conserved to the 
extent that it obeys the selection rule, g(W) = e 4iw ^ N W, corresponding to an R-charge of 2. 
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Thus, the invariance requirement of an order M superpotential monomial, W = Y\ I=1 $7, 
can be expressed by the condition, 

M M 

l(W) = KII = E $(1>i) + M = 2, (A.5) 

7=1 7=1 

implying the selection rule, J2f=i g(ipi) = 2 — M. 

Applying the above discussion to the S field dependent couplings, one derives the fol- 
lowing selection rules, valid for the three symmetry cases: 35 — 2 = 0, n + 2 + S = 
0; n + 2 7^ 0, 5 — 2^0, 2S — 2 ^ 0. It is again useful to single out the R like Peccei- 
Quinn symmetry, g PQ = P(R 2 A) n R 2 . This conserves all couplings with the exception of 
those involving the pair of Higgs boson superfields, for which one has the selection rules, 
~g n PQ (H d H u ) = 2 + n, ~g n PQ (H d H u S) = ~g n PQ (LH u S) = n + 1 + S. Focusing, for definite- 
ness, on the GLP, one obtains the following selection rules for the bilinear and trilinear 
couplings: g{H d H u ) = n + 2 ^ 0, g(LH u ) = m- n- p^0, g[LLE c ) = g(LQD c ) = 
m-2n-p-3 = -1, g(U c D c D c ) = m-2n-3 ^ -1, g(H d H u S) = S+n+l = -1, g(LH u S) = 
S + m-n-p-l = -l, g(S) =5-1^1, g(S 2 ) = 2S - 2 ^ 0, g(S 3 ) = 3S - 3 = -1. 

We can summarize the defining conditions for the GBP, GLP and GMP generators and 
the resulting representations of the generators by the formulas, 

GBP : m-2n-p-2 = 0; m-n-p + Sj^O, m - 2n - 2 ^ =>- g GBP = g PQ (RL) p 
GLP : m-2n-2 = 0, m-n-p + S = 0; m-2n-p-2^0 =>• g GLP = g PQ (L) p , 
GMP : m-2n-p-2^0, m - 2n - 2 ^ 0,m - n - p + S ^ =>• g GM p = ~g PQ R m ~ 2n ~ 2 U\ 

(A.6) 

which must be complemented by the conditions p ^ and n + 2 + S = 0. The anomaly 
cancellation conditions are now readily evaluated by inspection of the table given at the 
beginning of Section A which displays the fermion modes charges. One must include in 
the mixed gauge anomalies the contributions from the spin 1/2 gauginos, and in the grav- 
itational anomaly those from the gauginos and the spin 3/2 gravitinos which add to the 
anomaly coefficient 1 and —21 per mode, respectively. The contribution from the SM gauge 
group gauginos amounts then to J2 a dim(G a ) = 12. The anomaly coefficients for the gauge, 
gravitational and chiral anomalies are given by the formulas 

^ 3 = A{SU(3) 2 xZ N ) = 6- N g {4 + n)= rN, 
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A 2 = A{SU{2) 2 xZ N )=4- N g {A + n + p) + N 2h {n + 2) = rN, 
Al = A(U(lf Y x Z N ) = N 9 (-f - y + f ) + N 2h q + 1) = rN, 
A Z 2 = A(U(1) Y x Z 2 N ) 

= N g [l- 2(1 + mf + (1 - m + nf - (1 + n + p) 2 + (-1 + m + p) 2 ] 
+iV 2ft [(l - m + nf - (1 + m) 2 ] = riV, 



^aro* = A{grav 2 x Zjv) = iV fl (-15 - 5n + m - p) + A^ 2?i (4 + 2n) 

iV 

+12 - 21 + (-1 + S) = rN + ns—, 

2 

A z s = A(Z 3 N ) = N g [-Q - 3(1 + mf - 3(1 + n - mf - 2(1 + n + pf + (-1 + m + p) 3 ] 

iV 3 

+2AT 2h [(l + n - m) 3 + (1 + m) 3 ] + (-1 + Sf = rN + ns—. (A.7) 

o 

The selection rule, ( — 1 + S) = (— 1 + mx + + pz) = — (n + 3), may be used to remove 
the explicit dependence of the anomaly coefficients on the singlet field charges. The case of 
anomalous GS symmetries is analyzed in the same way as for the ordinary symmetries by 
introducing the shifted anomaly coefficients. The search of generalized parity solutions can 
also follow a similar two-stage procedure as described earlier. One solves in a first stage the 
anomaly cancellation equations at fixed iV for the integers (m,n,p), and in a second stage 
the selection rules n + 2 + S = 0, 3S - 2 = 0, 5-2^0, 2S - 2 ^ for the S field charges 
(x,y,z). To conclude, we note that the approach discussed here appears more systematic 
than the alternative one where one solves the anomaly cancellation equations after assigning 
charges to the various particles oq, au^, • • •, subject to the selection rules. 
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